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We interpret the forward Maxwell equation with up to third order induced polarizations 
and get so called nonlinear wave equation in frequency domain (NWEF) , which is based on 
Maxwell wave equation and using slowly varying spectral amplitude approximation. The 
NWEF is generalized in concept as it directly describes the electric field dynamics rather 
than the envelope dynamics and because it concludes most current-interested nonlinear 
processes such as three-wave mixing, four-wave-mixing and material Raman effects. We 
give two sets of NWEF, one is a 1+1D equation describing the (approximated) planar wave 
propagation in nonlinear bulk material and the other corresponds to the propagation in a 
waveguide structure. © 2013 Optical Society of America 

OCIS codes: 190.5530, 190.7110. 
Introduction 

There are three parts in the derivation. First, we review the derivation of Maxwell's wave equation in 
frequency domain and the expression of nonlinear induced polarization. Second, in the approximation of 
planar wave propagation, 1+1D nonlinear wave equation in frequency domain (NWEF) is derived as an 
interpretation of the 1+1D forward Maxwell equation with up to third order induced polarizations. Third, 
considering a waveguide structure, NWEF dependent on spatial mode profile is derived. 

1. Maxwell's wave equation and nonlinear induced polarization 

We start from Maxwell's equations and material equations shown below [IH5] : 

V* E = -f,VxH = -f +J ,1) 
V»B = 0, V«D = /9,V« J = (2) 



D = £ E + P,B = ^ H + M,J = (tE (3) 

"V" is the Laplace operator, E and H indicate the electric field (unit: ^) and magnetic field (unit: — ) 
vectors, D and B indicate the electric and magnetic flux densities. J is the current density vector and p 
is the charge density, representing the sources for the electromagnetic field. e$ is the vacuum permittivity 
(unit: ^) and /j,q is the vacuum permeability (unit: — ). P and M indicate the induced electric and magnetic 
polarization. 

In the absence of free charges in a nonmagnetic medium, we have M = 0, J = and p — 0. Hence the 
Maxwell's equations is simplified to: 
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<9B <9D 

VxE=--,VxH = - (4) 

After Fourier transforming the above equations (F(ui) = J_ F(t)e~ tut dt , matching the DFT algorithm 
in @MATLAB), we have: 

V x E = -iwB = -iufioH, V x H = iujt) = iw(e E + P) (5) 

To derive the wave equation regarding the electric field, we use the relation: 

V x (V x E) = V(V • E) — V 2 E (6) 

in which the left side can be extended as: 

V x (V x E) = V x (-ic^oH) = -iufm • (V x H) 

= -iujuq ■ iw(e E + P) = w 2 /i £oE + w 2 /i P = &oME + w 2 // P (7) 

On the right side, we have: 

V(V»E) - V 2 E« -V 2 E (8) 

which is supported if: 1) the high order induced electric polarizations (nonlinear induced polarizations) 
are considered perturbations to the first order induced polarization (linear induced polarization), i.e. P = 
Pl + Pjvl ~ Pl = 27T£oX ( ' 1 ' ) E; 2) the relative permittivity e r = 1 + 2t:x^ is independent on the spatial 
distribution, i.e. V • E = V • D/e e r = 0. 

Therefore, in frequency domain, we have the Maxwell's wave equation for the electric field: 

V 2 E + ^(w)E + w 2 ^ P = (9) 

In time-domain, it is: 

The induced polarization is further expressed as: 

p = p(l) + p(2) + p(3) + ... + p(m) ( U ) 

where the generalized expression of both the linear induced polarization P*- 1 ' = P^ and the nonlinear induced 
polarization p( m ) is: 

/oo />oo poo 

dh dt 2 --- dt mX (m Hti,t 2 ,--- ,t m )\E(t-t 1 )E(t~t 2 )---E(t-t m ) 
-oo J — oo J —oo 

/oo poo pOO 

dh dta-.- / dt mX {m} {t-t 1 ,t-t 2 ,--- ,t-t m )\E(t 1 )B(t 2 )---'E(t m ) 
-oo J —oo J —oo 

/OO pOO pOO 

duj! du2--- d Wm x (m) (wi,^2,--- ,w m )\E(w 1 )Efa)---E(u m )e it 'Z Ui (12) 
-oo J —oo J — oo 



= £o 



where: 



oo 



X^ m) (u> 1 ,u; 2 ,---u; m ) = — / dt x / dt 2 ---\ dfc m x (m) (*i > *a,---*m)c _iE,J ' t< (13) 



— oo J — oo 



is the temporal response function of the material, also called susceptibility in frequency domain, 
which is a (m+l)-rank tensor. The calculations among the electric fields arc dyadic product which result in 
an m-rank tensor. "|" indicate the multiple tensor product between two tensors, i.e. a (m+l)-rank tensor 
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and a m-rank tensor (or dyadic tensor). Therefore, the induced polarization p( m ) is a vector. In frequency 
domain, it has: 

/+oo 
p(m) e -iu,t dt 
-OO 

/OO pOO Z' + OO 

dwi ••• / dw m X (m) (wi, • • • ,w m )|E(wi) • • -E(w m ) / e-*^-^ 6 ")** 
-oo J—oo J — oo 

/oo poo 
doJi ■■■ / dw m x M (wi, • • • ,w m )|E(wi) • ••E(cj m )5(w - ^uj t ) (14) 
-OO J —OG 

Here, the delta function implies that the induced polarization always corresponds to the frequency which 
equals to the sum of the frequencies of the contributing electric fields. 

If setting r2 — i- e - ^ = d^i an d f^°° <K W — = 1, the above equation becomes: 

/OO /"OO /* oo 

dui / dw 2 • • • / dw m _i 
-OO «/ — OO «/ — OO 

, , x — ,m-l ^ — -m— 1 

x^'K^v',!!-^ ^)|E(w 1 )E( W2 )-.-E(fi-^ i (Ji),m>2 (15) 

and pW(fi) = 27re x (1) (^) »E(fi). 

Since p( m ) is a vector, it can be written as a sum of its components, each casting to one dimension, i.e. 
p(m) _ J^jPj™ 1 ^- Analogously, The nonlinear response tensor x^™"* has x^" 1 -* = Ei^-j i • Now is an 

,(m) 



-rank tensor and has RA = ^ 



(m) 



, where j, ai , • • • , a m is dimension mark. 



Then, the component of the induced polarization Pj is: 



/OO />00 />OG 

dfcj. / dt 2 --- dt m Rj TB) (t-ti,t-t 2j "-,*-tm)|E(ti)E(t 3 )"-E(t m ) 
-OO */ — OO */ — OO 

{/•OO />QO ^ 

/ tftl-"/ ^mx£;... am (*-*!,-•■ ,i-*m)-£ Ql (£l)---£ Qm (£m)l (1 6 ) 
J — oo J — oo J 



= £o 

Qi ■ ■ -a 

In frequency domain, it is: 



{/>oo /«oo m — 1 

/ duji--- I rfw m _ixj™ ) 1 ... Qm (wi, • • • ,<J m _i,fi - V Wi) 
J — oo J — oo 

x£; Ql (w 1 )---^ Qm _ 1 (w„ l _i)S Qm (r!-2j i Wi )Lm>2 (17) 

and P^. (Q) = 27T£ X) I Xj 1 ^ (fi)S 0l (0) >. xfcj...^ corresponds to one component of the tensor y^ m \ m 

which Xj-Lx * s coming from the matrix yp^ . 

Therefore, each component of the electric field Ej has a wave equation which, in frequency domain, is 
written as: 

V 2 4- + k 2 {oj)Ej + lo 2 ^P 3 = 0; P j =Pj 1] + Pj 2) + ■■■ + Pf n) (18) 
In particular, in uniaxial and biaxial crystals as well as cubic/isotropic materials, matrix x} 1 only has 
diagonal elements and therefore P^(f2) = 2ireox < j 1 j (^)Ej (ft) . we combine Pj with k^Ej and update the 
wave equation as: 

V 2 E 3 + k 2 (cu)E 3 + w*HoP jtI r L = 0; P jtNL = Pf + • • • + P { / n} (19) 
where k 2 (uj) = k'^l + 2nXj 1 j) and the refractive index in the dimension j is therefore defined as 
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1+1D NWEF 



In the approximation of planer wave propagation, we neglect the spacial dynamics in the propagation of the 
electric field but focus on the temporal and spectral dynamics, the Laplace operator is therefore reduced to 
the only derivative with respect to propagation axis z since the spacial dynamics is eliminated, i.e. V 2 — > ^-j. 
Therefore, Eq.(fT9|) is reduced to a 1+1D wave equation: 



d 2 ~ 







(20) 



By factoring out the fast dependence of the propagation coordinate from the electric field for all the 
frequencies, i.e. Ej(z,uj) = Aj{z,uj)e~ %k i^> z , we get: 



d 2 ~ d ~ 

g~2^j ~ 2ikj(w)—A 3 +uj 2 ^ P j ^Le tk ' z 







In the slowly varying spectral amplitude approximation (SVSAA) [4], i.e. 

-§^iAj <C -g^kjAj, which means the second-order derivative can be removed. 
So, we have: 



kjAj 



(21) 
we have 



dz Aj 



-i — , , Pj,NLe 



ikj z 



(22) 



and the reduced 1+1D Maxwell equation (also called forward Maxwell equation) regarding the electric 
field [5H7]: 



dE* 



ikj{uj)Ej 



-l — , - "j,NL 



(23) 



dz Jy ' J 2kj{u>) 
Among all the nonlinear induced polarizations, second-order and third order nonlinear induced po- 
larizations are most concerned. For second-order nonlinear induced polarization, the response is al- 
ways assumed instantaneous, giving rise to three- wave- mixing (TWM). The response function has 



Xj-.^aaC* 1 '* 2 ) = X^aiaa^C* 1 )^^)) where Xj^L ia2 is a constant indicating the response intensity. Correspon- 
dency, in frequency domain, the susceptibility has x^ aia2 (cJi, 0J2) = j^p^X < ^L 1 a 2 according to Eq.([l3]), i.e. 
constant for all the frequencies u>i and uj2- Hence, the second-order nonlinear induced polarization P- has 
the same form as shown in [51IH]: 



-(2) 



P?\t)=eoY,{xf aiQ2 E ai E a2 } 



(24) 



3V) 



so 



y (x (2) 



E 



}= £ ° E {xfl ia F[E ai E a A) 



(25) 



As for third-order nonlinear induced polarization, the response is not fully instantaneous but consists 
of the instantaneous electronic Kerr response and a fraction of non-instantaneous vibrational Raman re- 

R(h)6(t 2 - t!)S(t 3 ), 



sponse. The response function is therefore written as xj-onaaas^ 1 '* 2 '* 3 ' = ^ 
where R(t) = (1 — fn)5(t) + fnhn(t) and hn(t)dt = 1. fn indicates the amount of Raman frac- 
tion. hn{t) is the temporal Raman response function. The remaining instantaneous response is called 
electronic Kerr response which is the origin of the effects of self-phase modulation (SPM), cross- 
phase modulation (XPM) and four- wave-mixing (FWM). In frequency domain, the susceptibility is 



,(3) 



Y (3) 



(wi, w 2 ,w 3 ) 

(3) 



(2*-) 



W 3) 

6 Ai;q, 



(1 - /«) + /ft/ift(^i + W2) 



Hence, the third-order nonlinear induced 



polarization P- has: 
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P i 3) W = *> E {x5?aia a a, K 1 ~ h)E ai E a2 E a3 + f R (h R ® (i^I^))^] } 



(26) 



= e ° E {41^03 [(1 - fn)F [E ai E a2 E a3 ] + f R F [e^F' 1 [h R F [E ai E a2 ]]]]} (27) 



Finally, by substituting the nonlinear induced polarizations in Eq. (|23j) with their interpretations Eq.([25 
and Eq.(j27l), we get the 1+1D NWEF: 



dE 

—^-+ik j (ui)E j 
oz 



2c 2 k j (oj) 



E (4;«ia 2 f [E ai E a2 ] 



2c 2 kj(u 



E {x§Uo» [(! - /«) ^ [^Ba^a,] + fsF [e^F^ [h R F [E ai E a ,]]]]} (28) 



Equation (|28| can numerically be solved by split-step Fourier method together with Runge-Kutta method. 
It is noted that the frequency domain has a range (— oo, +00) and the contents of kj(cj) in negative frequencies 
are required. According to the causality of all the induced polarizations, kj(uj) shows a property of complex 
conjugate, i.e. its contents in negative frequencies are linked to what in positive frequencies. 



3. NWEF in waveguide structure 

In a waveguide structure, the spacial dynamics of the electric field is always dominated by the eigen modes 
in the waveguide, which determine the spacial distribution as well as the propagation constant of the elec- 
tric field. Since the orthogonality between any of the two eigen modes, we redefine the dimension mark 
j in the wave equation Eq. (fT9| to be the mode mark. Moreover, the electric field is redefined as [10] 
Ej(x,y, z,u) = Bj(x,y,uj)A^(z,uj) = Bj(x,y,uj)Aj(z,uj)e~' l ^ j ^ z , where Bj is the eigen mode distribution 
and fij indicates the mode propagation constant. Now Eq. (|19j) can be expanded as: 



\dx 2+ dy 2 + 1 Pj)^+^ e [ dz 2 lZ Pj dz l 



Aj + ui 2 jioPj^NL = (29) 



Remember all the eigen modes of the waveguide have (J^? + Jp- + k 2 — ft 2 ) = and using SVSAA, we 
get the reduced Maxwell equation for waveguide structure: 



f)F 



2j8i(w) 



(30) 



Compared with the 1+1D forward Maxwell equation Eq. (f2"3"|) . the only difference in the above expression 
is the replacement of the spacial propagation constant kj by mode propagation constant fij. However, it 
should be noticed that not only the dispersion characteristics but the nonlinear induced polarizations are all 
revised by employing the waveguide, which will be revealed later. 

Instead of focusing on the space-related electric field, now we concentrate on the spectral related amplitude 
Aj and get its dynamics by making spacial integral on both sides of Eq. (|3H)) . i.e.: 



1 2/3, 



dxdyB*P j:NL 



(31) 
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dA* 



+ i p j (u)A p j =-i— 



/ / dxdyB* Pj : NL, / / dxdy 

J J oo J J oo 



B, 



(32) 



Analogously, we should interpret the nonlinear induced polarizations in Eq. (|32|l in which the nonlinear 
susceptibilities are now space-dependent. For second-order nonlinear induced polarization, we employ Eq.(|17p 
and expand it with the definition of the electric field, getting: 



oo 

J J dxdyB*Pj 2 } = 27reo Yl J {duiAZ^A^iil-ui) 



x // dxdyx i ^l ia2 (x,y,uj 1 ,n-u 1 )B*(x,y,n)B ai (x,y,u 1 )B a2 (x,y,n~uj 1 )} (33) 

J J oo 

in which we can define a spacial integral factor as Phillips ct al. did in [10] : 



^QttaJ^ 1 '^ 2 ) = // dxdyx). i ! Xia2 (x,y,u} 1 ,uj 2 )Bj(x,y,(J-i + w 2 )B Ql (x, y, uJi)B a2 (x, y, u> 2 ) 

i OO 



1 



(2) 

dxdyx j ..a lQ2 (x,y)B*(x 7 y,uj 1 + ^2)B ai (x,y,u 1 )B a2 (x,y,u} 2 ) (34) 

(Z7T) JJoo 

and simplify the expression of second-order nonlinear induced polarization as: 

oo 

/ f dxdyB* = 2neo ]T [ du^^K^ ~ ^M^laM^ ~ ^) ( 35 ) 

JJco 

Now the spacial integral factor ®j-a ia2 pl a y s * ne role as the second-order susceptibility, determining the 
intensity of the nonlinear induced polarization. Moreover, it can be assumed that the variation of the mode 
distribution Bj with respect to frequency is slow enough compared with that of , or named slowly varying 
mode distribution approximation (SVMDA), therefore the spacial integral factor is considered constant 
©j 2 Q lQ , (i^i, ^2) ~ (2^®j 2 aia, an d the second-order induced polarization is further simplified to: 



dxdyB* Pf = £ 



Analogously, for third-order nonlinear induced polarization, we have: 



dxdyB* = e £ {eg^c, [(1 ~ fn)F [A^AIA] + f* F 



AfF- 1 



The spacial integral factor playing a role as the third-order susceptibility is: 



(36) 



(37) 



0j> 1 a 2 a 3 ( W l' W 2,W 3 ) 



Jj dxdyXj% ia20l3 (x, y, u)i, w 2 , u 3 )B* (x, y, u n )B ai (x, y, wi)B Q2 (x, y, w 2 )B a3 (x, y, w 3 ) 



1 - fn + fnhR(uJi + lo 2 ) 



(2tt) j 



// dxdyXj?l l0l2a3 Bj(x,y,y2 u n )B ai (x,y,Ui)B a2 (x,y,u) 2 )B a3 (x,y,u) 3 ) 
J J 00 ' " 

1-/b + /b^k(wi+W2)1 (38) 



(27T) 



-e (3) 

3 "j;aia 2 a 3 



ti 



Finally, by substituting the nonlinear induced polarizations in Eq. (|32[) with their interpretations Eq. (|36[) 
and Eq. ([37|) , we get the NWEF in the waveguide structure: 



- w,.,,, E {C,«..[< i -/*)- f, K- 4 !;,<J+/» F Kf- i [' i '! F K- 4 ™]]]]} < 39 > 
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